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Recent measurements of the mass of the top quark (t) are presented using 162 pb -1 of data of pp collisions at 
y/s — 1.96 TeV collected by the CDF detector at the Tevatron collider during Run II. The analyses focus on the 
semi-leptonic decay mode ti — » (blv)(bqq) with one or two identified bottom quarks (b). The Template Method 
reconstructs the invariant mass of the top quark in each event. The Multivariate Template Method enhances 
this approach by adding information on the event topology. The Dynamical Likelihood Method discriminates 
between possible mass values using top quark decay observables and attempts to use the maximum amount of 
information on top quarks provided by the Standard Model. All three methods produce similar results. The 
Dynamical Likelihood Method yields a top quark mass of 177.8±g;o (stat) ± 6.2 (sys) GeV/c 2 . 



1. INTRODUCTION 

The first experimental observation of the top 
quark (t) by the CDF and D0 collaborations at 
the Tevatron about ten years ago pp represented 
an impressive confirmation of the Standard Model 
(SM) and marked the beginning of a success- 
ful top quark physics program at the Tevatron. 
Properties of the top quark were first investigated 
using 110pb _1 of data at ^fs — 1.8 TeV collected 
in the Run I period between 1992 and 1996. The 
top quark mass, mt, was measured in all decay 
topologies {e.g. |2I3| ) arising from pp — > tiX fi- 
nal states. After a recent reanalysis of D0 Run I 
data 0] , the collaborations updated the combined 
Run I value to m t = 178.0 ± 4.3 GeV/c 2 0. 

The precise knowledge of m t is of great interest 
for particle physics. It is a fundamental param- 
eter of the SM, but more importantly, the ex- 
tremely high mass gives the top quark particular 
relevance in the calculation of other SM parame- 
ters. Electroweak corrections to the W propaga- 
tor introduce a quadratic dependence on the W 
boson mass mw and a logarithmic dependence on 
the mass of the long-hypothesized but still unob- 
served Higgs boson, m# . With the current preci- 
sion level of 2.5% for m t and < 0.1% for mw, the 
most likely value for mn obtained from global SM 
fits is m H = 114±^GeV/c 2 , and the upper 95% 
CL limit is 260GeV/c 2 jSj- A more accurate de- 
termination of mt would provide the opportunity 



for stringent consistency tests of the SM. 

The precision measurement of the top quark 
mass is therefore a very active topic in the current 
physics program of CDF in Run II. The upgraded 
Tevatron complex started in 2001, producing pp 
collisions at y/s~1.96TeV with steadily increas- 
ing instantaneous luminosities up to a record of 
~ 10 32 cm _2 s _1 (as of the time of this report). 
The higher collision energy led to an increase of 
the SM production ti cross-section of 35% [§]. 
The CDF detector j7j underwent substantial im- 
provements in the geometrical acceptance of the 
tracking and calorimeter systems as well as of 
the silicon vertex detector which is crucial for the 
identification of ILhadrons from t decays. 

All ti decay modes are being investigated in- 
tensively. Here we report on three recent CDF 
measurements of mt in the lepton plus jets chan- 
nel using 162pb _1 data collected between March 
2002 and September 2003. 

2. EVENT SELECTION 

According to the Standard Model, top quarks 
at the Tevatron are mainly produced in pairs 
through quark-anti-quark annihilation (85%) or 
gluon-gluon fusion (15%), followed by a prompt 
decay t — > Wb with a branching ratio BR^ 100%. 
The subsequent W decays define the event signa- 
ture. In the decay mode studied here, one W de- 
cays into two quarks, and the other W decays into 



an electron (e) or a muon (/i) and the correspond- 
ing neutrino (y). This channel has a manageable 
signal to background ratio S/B ~ 0(1) and a 
satisfactory BR ~ ^ that is viewed as a good 
compromise between all available decay modes. 1 

The final state has a spherical topology due 
to ti production near threshold and is character- 
ized by four jets, one charged lepton with high 
transverse energy (Et), and a sizable amount of 
missing transverse energy ($t) due to the unde- 
tectable v. Top quark candidate events are se- 
lected by requiring exactly one isolated e or /i 
with E T > 20 GeV and E T > 20 GeV. Further- 
more, at least three jets with Et > 15 GeV and 
one jet with Et > 8 GeV have to be present. The 
DLM analysis (explained below) requires exactly 
four jets with Et > 15 GeV. The jets must lie in 
the pseudo-rapidity region \rj\ < 2 so that they are 
well measured by the detector. Finally, one or two 
displaced secondary vertices (originating from the 
decay of long-lived B-hadrons) are identified by 
the silicon vertex detector ("6-tags"). This re- 
quirement improves the sample purity by a factor 
of 3-4 compared to the non-&-tagged case. 

Dominant background sources are Wbb and 
Wcc final states, non-W QCD processes and 
events where one jet is misidentified as a 6-jet. 
Minor sources of background are WW jW Z, Wc 
and single top quark (tb) events. 

3. TOP MASS MEASUREMENTS 

3.1. Template Method 

The Template Method (TM) closely follows the 
traditional mass measurement strategy in Run I. 
To get the momenta of the daughter particles of 
the assumed initial ti configuration, the raw jet 
energies are corrected for detector effects and ar- 
tifacts due to the limited jet cone size, as well 
as for hadronization and for various non-ii con- 
tributions of the same event. One part of the 
corrections are directly derived from the data, 
another part is derived using a Herwig ti Monte 
Carlo @ and a simulation of the CDF detector. 
The analysis treats the four leading Et jets as 
quark jets from the ti system. The event re- 

lr The dilepton mode (with e or fi) has S/B ~ 10, BR ~ ^j- , 
the all-hadronic mode has S/B < j^, BR ~ £. 
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Figure 1. Reconstructed top quark mass data com- 
pared with the combined signal plus background fit. 

construction is complicated by the combinato- 
rial problem of correctly assigning these jets to 
the four initial partons, yielding six (two) solu- 
tions ("permutations") for the single (double) b- 
tag case. The number of solutions doubles be- 
cause of the twofold ambiguity of the longitudinal 
component of the v momentum. 

For each configuration, an invariant top quark 
mass is reconstructed by minimizing a x 2 expres- 
sion calculated using the transverse momenta of 
the jets, the charged lepton, and the missing en- 
ergy. The x 2 takes into account the detector res- 
olution as well as the decay widths of the t and 
the W, and also incorporates various kinematic 
constraints (e.g. m t — mi, M w + = M w -). The 
solution with the smallest x 2 value is taken as the 
"reconstructed top mass" of the candidate event. 
This quantity is used to build templates for sig- 
nal Monte Carlo samples with varying generated 
mt but also for various background samples. The 
signal templates are parametrized by continuous 
functions of m t . The templates serve as prob- 
ability densities to be compared with the data 
in order to derive a top quark mass probability 
for each event. The most likely value for mt is 
extracted using an unbinned likelihood fit that 
maximizes the joint probability density (= prod- 
uct of all event-by-event probability densities) for 
the whole sample of reconstructed masses. It also 
includes a background fraction constraint. 

With a sample of 28 candidate events (Fig. ^) 
and 6.8 ± 1.2 expected background events 
the method yields m t — 174.9^y'Y(stat) ± 
6.5(sys) GeV/c 2 . The systematic error is dom- 
inated by uncertainties of the jet energy scale 
(±6.3GeV/c 2 ) which mainly arises from the cal- 
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ibration of the CDF calorimeter to a uniform re- 
sponse and setting the absolute jet energy scale. 

3.2. Multivariate Template Method 

The Multivariate Template Method (MTM) 
has three main features: 1) The jet energy scale 
(JES) is adjusted event-by-event. 2) Three types 
of signal templates are used to handle the combi- 
natorial problem. 3) Other kinematic variables, 
in addition to mt, are used to improve S/B dis- 
crimination. The method performs the same 
generic jet corrections as TM but uses a refined 
kinematic fitter for the event reconstruction. In 
particular it includes a JES factor calibrated by 
the reconstruction of the W — > qq' decay using 
a W mass constraint. On average this partially 
compensates for systematic JES shifts but also 
increases the statistical error because of event-by- 
event fluctuations of the JES. The fitter contains 
a JES constraint which allows the balance of JES 
systematics vs. statistical uncertainties. Further 
kinematic constraints similar to TM are imposed. 
Only the neutrino solution for which the t and t 
masses are closest is taken. 

Three signal templates are used to handle jet- 
parton combinatorics: correct permutation sam- 
ples, incorrect permutation samples, and incor- 
rect jet assignment samples (where e.g. a gluon jet 
is assigned to a quark). It can be shown that the 
knowledge of the event template type improves 
the mass resolution by a factor of ~ 1.7 in an ideal 
scenario with no background. MTM attempts to 
predict the correct permutation probability by us- 
ing the differences between the smallest x 2 an( i 
the x 2 values of the other permutations, and by 
using kinematic information with weak m t depen- 
dence, specifically the angle between the lepton 
and 6-quark in the W — > Iv rest frame, and tt 
spin correlation variables. These probabilities are 
used to weight the signal templates accordingly. 

In addition to the reconstructed top quark 
mass, further kinematic variables with high S/B 
separation power are used to construct multidi- 
mensional templates. From studies of various 
histogram divergence measures it turns out that 
the scalar sum of the transverse momenta of the 
four leading jets (Pt4) is a promising variable to 
distinguish signal from noise. The templates are 
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Figure 2. Reconstructed mass vs. Pt4 for the data 
and for signal and background expectation. 

generated using Kernel Density Estimation JHj , a 
non-parametric density reconstruction technique 
particularly suited for multivariate densities. 

The top quark mass likelihood is expressed 
in terms of multivariate probability densities for 
signal and background and considers the three 
types of signal templates. The background frac- 
tion is allowed to float freely. The JES con- 
straint is optimized w.r.t. the expected total er- 
ror. A maximum likelihood method applied to 
a sample of 33 selected events (Fig. |2J yields 
m t = 179.6l|;f(stat) ± 6.8(sys) GeV/c 2 and a 
most probable background fraction of 0.34 ±0.14. 

3.3. Dynamical Likelihood Method 

The Dynamical Likelihood Method (DLM) is 
an original CDF method JT] that attempts to 
use the full amount of information about the tt 
process provided by the Standard Model. The top 
quark mass likelihood L^(m t ) for the i th event 
is given by 

L«(m t ) = lY. E^l^l^'^/bt) 

J * — ' * — ' flux 
perm zysol. 

xio(x,y;mt) dx 

where x denotes the parton momenta and y de- 
notes jet observables. F{z\, zi) is the parton dis- 
tribution function for the incoming partons, and 
f(pt) is the probability that the tt system has 
received a transverse momentum p t due to ini- 
tial state radiation. The matrix element M for 
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Figure 3. DLM joint negative log likelihood. 

the tt production and decay process is linked to 
the four-jet configuration using transfer functions 
w(x, y; m t ) (TF). These are Bayesian probabili- 
ties that a tt parton configuration x with a given 
mt was generated when y was reconstructed. The 
TFs are functions of (E^ Iton - E^)/E^ ton and 
parametrized in jet Et and r\ bins. They are cal- 
culated for W or b jets separately using a Herwig 
tt MC for a fixed top quark mass. DLM performs 
some of the jet corrections as TM and MTM {e.g. 
relative calorimeter response, calorimeter to par- 
ticle correction), while it handles the transition 
from partons to particles via TFs. 

The likelihood function LW sums all per- 
mutations and neutrino solutions. The most 
likely top quark mass is obtained by maximiz- 
ing IljLW(mt). The joint negative log likeli- 
hood curve extracted from a sample of 22 se- 
lected events is shown in Fig. |3l It should be 
noted that the likelihood considers only the sig- 
nal but not the background matrix elements. The 
background contribution is hence minimized by 
requiring exactly four jets. Residual background 
effects, but also the top quark mass dependence 
of the TFs, are taken into account by applying 
an MC-based correction ( "mapping" function) to 
the maximum likelihood mass. Using 4.2 ± 0.2 
expected background events, this method yields 
m t = 177.8lg:o(stat) ± 6.2(sys) GeV/c 2 . 

4. SUMMARY AND OUTLOOK 

Three recent CDF measurements of the 
top quark mass have been performed in the 
ti lepton plus jets channel using 162pb _1 of 
Tevatron Run II data. The results are m t = 



174.9l^(stat) ± 6.5(sys) GeV/c 2 (TM), m t = 
179.6^;3 (stat) ± 6.8(sys) GeV/c 2 (MTM), and 
m t = 177.8t^(stat) ± 6.2(sys) GeV/c 2 (DLM). 2 
The systematic error is dominated by jet energy 
uncertainties. MTM allows for an optimal com- 
promise between statistical and systematic errors 
and will be more powerful at higher accumulated 
luminosity. The DLM result has a considerably 
less statistical uncertainty than the other meth- 
ods and is competitive with the best CDF Run I 
measurement 0. 

Tevatron is performing well and is expected to 
provide 400 pb" 1 by the end of 2004, lftT 1 by 
the end of 2005, and 4.4-8.5 fb" 1 by the end of 
Run II, thus leading to a significant improvement 
of the statistical error. A priority is the reduc- 
tion of the jet energy systematics. The stated 
CDF goal for the total error is 3GeV/c 2 . Sub- 
stantial progress is expected soon, e.g. due to re- 
cent developments in the calorimeter simulation. 
Furthermore there is a substantial potential to re- 
fine the analyses. Through these measurements, 
CDF has established techniques that will provide 
a precise determination of the top quark mass in 
the future. 
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2 The results are consistent with an updated CDF mea- 
surement in the dilepton channel reported after this con- 
ference: m t = 176.5+jgg(stat) ± 6.9(sys) GeV/c 2 . 



